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The public code GOSAM for the computation of the one loop virtual corrections to scattering 
amplitudes in the Standard Model and beyond is presented. Particular emphasis is devoted to the 
interface with other public tools via the Binoth Les Houches Accord. We show with examples 
that doing LHC phenomenology including automatically Next to Leading Order QCD corrections 
is now handy. 
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1. Introduction 

The recent discovery of a new resonance at the Large Hadron Collider (LHC) at CERN with 
mass 125 GeV [1] is the result of the efforts in the construction and operation of the accelerator 
and of the detectors as well as of the efficiency in the implementation and interpretation of very 
detailed analyses on the large amount of informations collected. Now the target is to perform de- 
tailed measurements of all the branching ratios of the novel particle while searching for further 
new resonances in all the accessible mass regions. The first question is to establish if the new res- 
onance is or not the Higgs boson. In this scenario the help of automated tools for the simulation 
of signals and backgrounds coming from the Standard Model or beyond is desirable and is now 
possible including Quantum Chromo Dynamic (QCD) corrections up to the Next to Leading Order 
(NLO) accuracy in the strong coupling constant. The QCD at Hadron Colliders play of course the 
main role and several fully automated tools exist that include the leading order matrix element, the 
leading logarithmic correction as a parton shower and the subsequent hadronization, leading to full 
event simulation [2]. The NLO hard QCD corrections, being in general quite sizable, also have to 
be included to get predictions with a theoretical error comparable with the accuracy attainable at 
the LHC. NLO QCD corrections have been successfully matched in full generality with a parton 
shower avoiding the double counting that derive by the fact that the shower produces also some of 
the emission predicted in the full NLO calculations [3] . In this way a number of processes have 
been studied at NLO accuracy plus parton shower. The deep understanding of fundamental prop- 
erties and the structure of the divergences of QCD corrections up to the NLO made the automation 
of the subtractions of the divergences and the matching to the parton shower feasible in full gener- 
ality. On the other hand progress has been done in the computation of the virtual part of the NLO 
correction exploiting the analytic properties of loop amplitudes. Advanced computer languages 
and the increasing cpu performances give now the possibility to automatize also the computation 
of the virtual corrections even for complicated hard processes, i.e. processes with many legs in the 
final state. 

A key ingredient for the full automation of NLO computation is given by the Binoth Les Houches 
Accord [4] (BLHA) that fixes a set of universal rules for the communication among different pro- 
grams that can provide the complementary ingredients needed for the full event simulation. 
It is the aim of this presentation to illustrate the usage of GOSAM [5], a new fully automated pack- 
age for the computation of the one loop matrix elements, and its interfacing with other programs for 
the integration over the phase space of both the real and the virtual corrections to Born processes. 

2. GOSAM 

GOSAM is a python package that has been designed for the generation of codes for the numer- 
ical evaluation of one loop amplitudes. At the moment the code produced are in Fortran95. Inside 
GOSAM a number of libraries are used: the QGRAF [6] Fortran library is used for the generation 
of the Feynman diagrams, while FORM [7] is used for the algebraic treatment of the diagrams 
mainly through the SPINNEY [8] FORM package, the formulas are then optimized and written by 
the java package HAGGIES [9], finally once a code for a process has been generated it runs calling 
libraries for the reduction of tensor integrals and scalar integral evaluation [10]. 
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The construction of the codes has a very simple structure, the aim is to write efficiently the numera- 
tor of each Feynman diagram (or a group of them) as a multivariate polynomial of the components 
of the loop momentum. For each physical phase space point such polynomial in the loop mo- 
mentum gets numerical coefficients related to the specific kinematic configuration and as a tensor 
integrand sitting on a certain set of Feynman denominators is reduced and computed numerically 
by dedicated libraries. 

After several years since the seminal work of Passarino and Veltman [11], loop tensor integral 
reduction techniques have seen new progress. The original algorithm has been extended to op- 
timize the numerical evaluation of the form factors especially near unstable configurations [12]. 
On the other hand, the parametric form of the residues of the multipole expansion of the Feyn- 
man integrands has been fully determined and exploited to construct recursive algorithms based on 
polynomial interpolation, leading to the direct numerical determination of the coefficients of the 
scalar integral [13]. 

In general two factors determine which method is the fastest or more efficient: the efficiency of 
the tensor reduction library and the way the numerator is prepared and matched with the tensor 
reduction library. Further, having more then one construction leads to valuable consistency check 
on the results. Alternative methods for the construction of the numerator function to match with 
loop integrals have also been proposed in [14]. A schematic representation of the functionality of 
the codes generated by GOSAM is given in Figure 1. The basic tensor reduction implemented is 



Figure 1: Reduction flow of the GOSAM codes for the evaluation of the numerical virtual amplitudes. 

through the SAMURAI library, based on the integrand reduction method mentioned above, that 
also provides tests to establish the quality of the result. If the test is not passed then the reduction 
of the tensor integrals is repeated with GOLEM95. For these cases the codes are equipped with a 
second tensorial representation for the numerator and give more accurate results while remaining 
in double precision. The reduction through SAMURAI is the default because it turns out to be a 
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Table 1: Sample of processes for which GOSAM has been compared to the literature. 



process 


process 


e + e~ — > uTi 


pp — > /'/' 


e + e~ — > tt 


W (massive b's) 


uu — ¥ dd 


e+e~ e+e~7(QED) 


88^88 


pp — > ttH 


88-^8^ 


pp tiZ 


pp^rtt 


yy -» 7777 (fermion loop) 


bg^Hb 


pp -> 


yy -> yy (W loop) 


— > £>M>£> 


pp^W ± ./(QCD corr.) 




pp->W ± ;(EWcorr.) 


pp — > ffM> 


pp -> W ± t 





bit faster then the one through GOLEM95. 

This switch happens dynamically during the computation of the virtual matrix element, on a dia- 
gram by diagram (or group of diagrams) base. Finally, for QCD corrections, the renormalised result 
for the single pole is compared to the expectation given by the well known formulas expressing the 
universal behaviour of the infrared QCD divergences, and if the agreement do not fulfil the required 
accuracy the whole amplitude is computed again using GOLEM95. We have not performed a full 
study of the correlation among the numerical instabilities of our procedures and the relative impor- 
tance of the virtual correction. Experimentally we find that fixing the threshold for the mismatch 
of the single pole as 

Single pole(Exact) — Single pole{Numerical) 4 
£= WBom ~ ° ( } 

has two consequences: first we note that only a very small fraction of points turns out to be above 
threshold (typically one per million for two to three processes) and, second, the results for the 
distributions turn out to be in excellent agreement with the one obtained with codes containing 
analytic formulas for the virtual amplitudes like the MCFM [15] suite of programs. Finally, if this 
higher level instability trigger finds a phase space point to be still unstable this can be saved on a 
file to be reprocessed with higher precision with the same GOSAM code that can be compiled even 
with intermediate or quadruple numerical precision. 

A number of pointwise comparisons has been performed with results obtained mainly with other 
codes or by coding amplitudes reported into the appendix of the relative papers, or simply compar- 
ing the numerical results reported for benchmark phase space points. A sample of the comparisons 
is given in Table 1 , while in Table 2 a list of timings is reported for the generation of the codes and 
the relative running time per phase space point. 

We generated codes for the NLO QCD corrections for the production at hadron colliders of four 
bottom quarks [16] and for the associated production of two W bosons with opposite charge plus 
two jets [17]. Both computations allowed to test the packages on amplitudes consisting of roughly 
one thousand diagrams. The running time per phase space point for such matrix elements is around 
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Process 


Generation [s] 


Evaluation [ms] 


bg^Hb 


236 


2.49 


dd it 


324 


4.05 


dg^rdg 


398 


3.08 


e + e- -> it (LanHEP) 


180 


1.27 


e + e~ — > mm (AutoTools) 


173 


0.64 


<?<? (LanHep) 


1022 


1.70 


gg^gZ 


529 


15.18 


gg^tt (UFO) 


1225 


29.45 


H -^yy 


140 


0.24 


gb^e-v e t 


337 


2.89 


ud — > e~ V e 


71 


0.09 


ud — > e~V e g 


154 


1.15 


uu —> dd 


186 


2.06 


ud -> W + W + cs 


1295 


17.37 


yy^yy 


597 


6.08 



Table 2: Time required for code generation and calculation of one phase-space point. The results were 
obtained with an Intel(R) Core(TM) i7 CPU 950 @ 3.07GHz. The time for the evaluation of a phase space 
point is taken as the average of the time obtained from the evaluation of 100 random points, where the code 
was compiled using gf ortran without any optimisation options. The generation of the R2 term was set to 

explicit. 

one second for a single subprocess preventing the direct integration on a single core of the full 
virtual matrix elements. In these cases we proceeded through the re-weighting of a previously 
obtained sample of Born level un-weighted events to obtain a reasonable estimation of the virtual 
correction. The integration of the real part and of the subtractions has been performed with the 
packages MadGraph [18], MadDipole [19] and MadEvent [20]. The gluon initiated four bottom 
quarks production virtual amplitudes were integrated for the first time using for their numerical 
evaluation a GOSAM code, while the main contributions for the process pp —> W + W~ + 2 jets 
were already be computed in [21]. In [17] we included the computation of the closed fermion 
loops with electroweak gauge bosons attached and the non top resonant contribution form the third 
generation. 

3. Interfacing GOSAM 

Pointwise tests and integrations based on the re-weighting procedure for complicated two to 
four one loop processes show how our public package GOSAM can be efficiently used as a fully 
automatic generator for the virtual amplitudes. However, that is only a single piece in the task 
of the full automatic NLO code generation. The standard communication setup established in the 
BLHA gives then the possibility to link GOSAM to any other program that supports the same 
standard and supplement the other ingredients. The level of the automation in the generation of 
the tree level Born and real corrections, of the subtractions of the divergences as well as in the 
management of the integrations and the analysis of the results is already very high. We exploited 
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the standard to link GOSAM to the SHERPA package that since the release 1 .4 is equipped with 
all the ingredients needed for the parton level next to leading order simulation of QCD correction 
to the Standard Model processes, so that the whole NLO code is generated by a single command- 
line executable [22]. In this way the generation of the whole code is fully automatic so that once 
GOSAM provides the virtual amplitudes all the features already implemented in SHERPA for the 
leading order analyses become available for the NLO computation. Representative examples of 
distributions obtained with code generated fully automatically with GOSAM+SHERPA are given 
in Figures 2, 3 and 4 for the processes pp — > W~bb with massive bottom quarks, pp — > W~ + 2 jets 
and pp — > W + W + + Ijets respectively. The comparisons have beed done with MCFM [15] for the 
former two processes and with [23] for the latter. In conclusion we have presented the GOSAM 
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Figure 2: Charged lepton and W transverse momentum distribution for the process op- 
tion obtained with GOSAM+SHERPA are compared to MCFM for a set of standard cuts. 



W bb. Predic- 




Figure 3: Leading jet transverse momentum and diet invariant mass distributions for the process pp — > 
W~ + 2 jets. Prediction obtained with GOSAM+SHERPA are compared to MCFM for a set of standard 
cuts. 

package for the generation of numerical codes for the virtual amplitudes at one loop and shown 
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W + W + + 2 jets: leading jet rapidity W + W + + 2 jets: second leading jet transverse momentum 
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Figure 4: Rapidity of the leading jet and second jet transverse momentum distributions for the process 
pp —> W + W + +2jets. Prediction obtained with GOSAM+SHERPA are compared with plots extracted 
from [23] following the phenomenological analysis reported there. 



how the interface with other public codes can be exploited to have the fully automatic generation 
of programs for the computation of the NLO correction to the hard scattering processes. 
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